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Abstract 



Supersymmetric (SUSY) Ward identities are considered for the N=l SU(2) 
SUSY Yang Mills theory discretized on the lattice with Wilson fermions (glu- 
inos). They are used in order to compute non-perturbatively a subtracted 
gluino mass and the mixing coefficient of the SUSY current. The compu- 
tations were performed at gauge coupling (5 — 2.2) and hopping parameter 
«; = 0.1925, 0.194, 0.1955 using the two-step multi-bosonic dynamical-fermion 
algorithm. Our results are consistent with a scenario where the Ward identi- 
ties are satisfied up to 0{a) effects. The vanishing of the gluino mass occurs 
at a value of the hopping parameter which is not fully consistent with the 
estimate based on the chiral phase transition. This suggests that, although 
SUSY restoration appears to occur close to the continuum hmit of the lattice 
theory, the results are still affected by significant systematic effects. 



1 Introduction 



A better understanding of non-perturbative phenomena in supersymmetric (SUSY) gauge 
theories could be gained in the framework of the lattice regularization. An immediate 
difficulty arises, however, because the lattice regularized theory is not supersymmetric 
as the Poincare invariance, a sector of the superalgebra, is lost. This is evident if one 
considers the super-algebra in the canonical formalism (the notation is with Weyl spinors) 

{Q,,Q4 = 2<^P, . (1) 

This relation cannot be fulfilled in a discrete space-time manifold, where momenta are 
not generators of infinitesimal space-time translations. More specific difficulties arise in 
the fermionic sector where spurious states may violate the balance between bosonic and 
fermionic degrees of freedom. In the standard approach to lattice gauge theories with 
Wilson fermions (for an approach using do main- wall fermions see [jl|) the suppression of 
spurious states in the continuum limit is obtained by adding to the action an 'irrelevant' 
term (Wilson term) which explicitly breaks SUSY. 

Some time ago Curci and Veneziano proposed that in spite of this substantial 
SUSY breaking the Wilson discretization may be safely applied to SUSY gauge theories: 
the symmetry is recovered in the continuum limit by properly tuning the bare parameters 
of the action. They considered the simple example of the N=l supersymmetric Yang Mills 
theory (SYM). This is the supersymmetrized version of quantum gluodynamics where the 
A*"^ — 1 gluons are accompanied by an equal number of fermionic partners (gluinos) in 
the same (adjoint) representation of the color group. The Wilson action for the N=l 
SYM breaks SUSY by the Wilson term for the gluino action and by a gluino Majorana 
mass term. The lattice Ward identities (Wis), considered in |^ for the case of the chiral 
symmetry in QCD (see also 0, 0), provide a general theoretical framework for properly 
dealing with the problem of the restoration of the symmetry in the continuum limit. As a 
consequence of the SUSY breaking, similarly to QCD, the gluino mass is shifted: a tuning 
procedure on the bare mass is required in order to recover massless gluino and SUSY in 
the continuum limit. 

The problem of the Monte Carlo simulation of the N=l SU(2) SYM with Wilson 
fermions was considered by this collaboration in the past p, |^, ||, |, |10|] (for a study in 
the quenched approximation see [0). The theory with dynamical gluino was simulated 
by the two-step multi-bosonic (TSMB) algorithm defined in |T^, |TB[. The formulation of 
the algorithm, flexible with respect to the parameter Nf, allows to treat the gluino which 
is a Majorana fermion (A^/ = l/2). An extensive analysis of the low-energy aspects of the 
SU(2) SYM was performed. 

An intricate point in the N=l SYM is the tuning of the theory to the massless gluino 
limit. The straightforward spectroscopic method familiar in QCD cannot be applied in the 
case of the N=l SYM since the theory possesses only an anomalous axial chiral symmetry 
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and no Goldstone boson is available. The problem was considered in by studying the 
behavior of the finite-volume gluino condensate as a function of the hopping parameter: 
the massless gluino is expected to correspond to two degenerate vacua with symmetric 
probability distribution of the gluino condensate. SUSY restoration can be also verified by 
direct inspection of the low-energy mass spectrum [^: this is expected to reproduce the 
SUSY multiplets predicted by the low-energy effective Lagrangians |T^, |T5[. An accurate 
analysis of the spectrum is however a non-trivial task from the computational point of 
view and an independent method for checking SUSY is welcome. A possibility |]TT| is to 
use the SUSY Wis to determine the gluino mass in the same way the chiral Wis are used in 
QCD to determine the quark mass. A great simplification consists in considering the on- 
shell regime. In addition the WI approach improves the insight in the renormalization of 
the lattice SUSY current. The properly renormalized current defines the supercharge and 
satisfies the appropriate superalgebra. The renormalization of the lattice SUSY current is 
more complicated compared to the chiral case since SUSY is more severely broken on the 
lattice. Explicit one-loop calculations in lattice perturbation theory may lead to a better 



understanding of this problem [16, 17 



In this work we concentrate on the SUSY WI approach. As we shall see, two unknown 
parameters appear in the on-shell SUSY Wis. These are ratios of three coefficients enter- 
ing the Wis: Zs and Zt multiplying the divergences of the SUSY current df^S^{x) and of 
the mixing current d^T^{x), and the subtracted gluino mass m^. For ms = 0, SUSY is 
expected to be restored in the continuum limit. As a result of the study we obtain a non- 
perturbative determination of the dimensionless ratios amsZ^^ and ZtZ^^. We consider 
the gauge coupling (3 = 2.3 on a 12^ x 24 lattice and three values of the hopping parameter 
K = 0.1925, 0.194, 0.1955, corresponding to decreasing gluino mass. Preliminary results 
were presented in ||18|, |19|, 

The paper is organized as follows. In Sec. 2 the Curci-Veneziano approach is intro- 
duced in the case of the N=l SYM. The lattice action is defined and its symmetries are 
pointed out (see also Appendix A). The latter play a special role in the analysis of the 
SUSY Wis. In Sec. 3 the formalism of the lattice Wis is reviewed. The fundamental issue 
in this context is the renormalization of the 'irrelevant' operator entering the Wis because 
of the explicit breaking of SUSY in the lattice model. A detailed account is given in Ap- 
pendix B with an analysis based on the discrete hypercubic group. The result is that the 
SUSY Wis assume a specific lattice form. We discuss in this paper only the simplified case 
of the on-shell regime. Suitable gluino-glue insertion operators are discussed in Sec. 4. In 
Appendix C the role of the symmetries in this context is clarified. In Sec. 5 we give an 
account of the present setup of the TSMB algorithm. The parameters have been tuned 
in order to get good performance for light fermionic degrees of freedom. We also measure 
some quantities (smallest eigenvalue, sign of the Pfaffian, the scale tq) which characterize 
the set of configurations under study. Sec. 6 is devoted to the numerical analysis of the 
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SUSY Wis with an outline of the method and the presentation of the numerical results. 
Conclusions are finally drawn in Sec. 7. 



2 The N=l SYM on the lattice 

We adopt the formulation of Ref. for the lattice discretization of the N=l SYM with 
Nc colors. The pure gauge action Sg is the standard plaquette one 



where the plaquette is defined as 



-T" I I— /-I I ^ ' 



(2) 



= Vlix) Ulix + 0) U,{x + (i) U^{x) , (3) 

and the bare gauge coupling is given hj (3 = 2Nc/gQ. For Wilson fermions the action 
(with Wilson parameter set to r = l) reads 

Sf = 5]a^Tr[l5^(A(x)(7,-l)f/;(x)A(x + /i)f/,(x) 

- X{x + fi) (7^ + 1) U,{x) \{x) Ul{x)) + (mo + ^) \{x)\{x) ] , (4) 

with a the lattice spacing and mg the gluino bare mass. The gluino field A(x) is a 
Majorana spinor transforming according to the adjoint representation of the gauge group. 
The symbol 'Tr' denotes the trace over the color indices. In this work we consider Nc = 2, 
for which the adjoint gluino field is expressed in terms of Pauli matrices ar 

3 1 
■r=l 

The following relation (Majorana condition) holds for an Euclidean Majorana field: 

A(x) = \^{x) = CA^(x) , (6) 

where C = 7072 is the spinorial matrix associated with the charge-conjugation symme- 
try C. Boundary conditions are taken to be periodic except for fermionic fields in the time 
direction, which are anti-periodic. 

In Monte Carlo simulations a different parametrization is used. The hopping param- 
eter K is defined as ^ 

^ 2(4 + moa) 

and the fermionic action is expressed in terms of the fermion matrix Q 

= f^Y.^\y)Qys,xry{x) (8) 



2k 

xr,ys 
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(Dirac indices are implicit). The fermion matrix is given by 

4 

Qys,xr — Q ys,xr\U^ = Syx^sr ^ ^ ^ ~l~ T/^)^r,2:/i ~l~ ^y+fi,x{^ 0'/i)^r,j//i] l'^) 

with the adjoint hnk Vrs^x^i{x) defined as 

Vrs,x^, = Vrs,x,,[U] = ^ Tr (f/^^d^t/^^cT,) = V;^^x^ = V-^J^ . (10) 

The action (@)-(§) is invariant under the discrete symmetries P (parity), T (time- 
reversal) and C (charge-conjugation). For the case under consideration the latter symme- 
try implies the following relation for the fermion matrix 

Ql,ys[U] = CQy,,xr[U]C-' (11) 

(transposition is intended on the suppressed Dirac indices). The discrete symmetries of 
the lattice action play an important role in the subsequent analysis of the Wis; their 
explicit definition is given in Appendix A. 



3 Lattice SUSY Ward identities 

SUSY is explicitly broken in the action (0)-(|5) by the gluino mass term, by the Wilson 
term and by the lattice discretization. Using lattice SUSY Wis a soft-breaking subtracted 
gluino mass ms can be defined. The expectation is that the vanishing of ms', for asymptot- 
ically small lattice spacings, ensures the restoration of SUSY up to discretization effects. 
In this Section we discuss these issues, which were first introduced in and have also 



been considered in WHA and 



Lattice SUSY transformations complying with gauge invariance, P, T and the Majo- 
rana nature of the gluino field are ||2l], |T6|Q- 

6Ui,{x) = -^^(^(x)7^[/^(x)A(x) + e{x + fihf.Xix + fi)U^{x)) , 

K(^) = '-^ie{xh,x{x)ul{x) + e{x + f,h,ul{x)x{x + i,)) , 

SX{x) = ^P^l^\x)a,J{x) , 

6X{x) = -le{x)a,.Pl^\x) , (12) 

where d{x), 6{x) are infinitesimal Majorana fermionic parameters. The lattice field tensor 
pIw\x) is clover- symmetrized so as to comply with P and T: 



^Our definition of the link variable Ufi{x) differs from that of (see our definition of the plaque- 
tte (0)); the two definitions are related by Hermitian conjugation. 



where 



U;^J{x) = Ul{x)Ul{x + u)U,{x + fi)U,{x) = U,,{x) , 

US}i^) = Ul{x)U,{x + A - 0)U,{x - u)Ulix - u) , 

Uj:n){^) = Uu{,x - u)Uy{x - fi- i>)Ul{x - fi- i>)Ul{x - jl) , 

UlS{x) = Upix - jj)Ul{x - ij)Ul{x -fi + u)U,{x) . (14) 

For any operator Q{y) the expectation value {Q{y)) is invariant if in the functional 
integral a change of variables according to the above SUSY transformations is performed. 
In the case of a gauge invariant operator Q{y) this results in the following WI 

E((V/^^r(^))^(^)> = rno {x{x)Qiy)) + {X^''\x)Q{y)) - (^^y (15) 
The SUSY current Sjl"^\x) is point-split {ps) |]16[ 

Sjf'\x) = -i5^a,.7MTr(p^(:')(a;)f/t(a;)A(x + /i)f/,(a;) +p(:')(a; + ^)t/,(x)A^^^ , 

pa 

(16) 

and the lattice derivative is the backward one V|^/(a;) = (/(x) — f{x — ft)) /a. We recall 
that SUSY is broken by the presence of a non-zero bare mass in the action, by the Wilson 
term and by the discretization. The first type of SUSY breaking gives rise to the term of 
the WI (|T^) involving the operator x{x)'- 

x{x) = Y,a,^TT{Pl,f{x)\{x)) . (17) 

pa 

The rest of the SUSY breaking results in the presence of the X'^^^\x) term. Its exact 
expression is not needed in the following. It suffices to know that in the naive 
continuum limit Xs{x) ~ aOii/2(a;), where Oii/2(a;) is a dimension- 11/2 operator. 

The last term in (^) is a contact term which vanishes when the distance \x — y\ 
is non-zero. This corresponds to the on-shell situation. We shall restrict ourselves to 
this regime in the numerical analysis of the Wis and contact terms will be consequently 
disregarded in the following discussions. 

The definition of the SUSY current on the lattice is arbitrary up to terms which vanish 
in the continuum limit. Another choice is the local (/oc) current 

Sjl-'^Hx) = -J2 ^p.7MTr (p(^')(^)A(a:)) . (18) 

pa 

This definition is preferable on the classical level |^ and is more convenient for analytic 
perturbative calculations. The local current Sp°^\x) satisfies a WI of the form (|^), with 
a symmetric lattice derivative V^/(x) = (/(x + fi) — f{x — fi))/2a (required to preserve 
P and T) and a SUSY-breaking term XC"'^) = X^f") + 0(a). 



3.1 Renormalization 



The WI (^) is a relation between bare correlation functions. The role of the symmetry- 
breaking operator X^p^\x) (or X^''°^\x)) is of particular interest, as it is related to current 
normalization and gluino mass subtraction. Its treatment in the present case follows 
closely that of the axial Wis in QCD 0, ^ . 

We consider the renormalization of the dimension- 11/2 operator On/2{x). According 
to the usual prescriptions, this implies mixing with operators of equal or lower dimensions 
d, which have the same transformation properties under the symmetries of the lattice 
action. A discussion of the mixing pattern on the basis of the discrete hypercubic group 
is carried out in Appendix B. The result is that no Lorentz-breaking mixing arises at least 
in the on-shell regime. The mixing pattern (involving operators with 7/2<d < 11/2) in 
the on-shell case is given by 

Ofi/2(x) = Zn/2[On/2(x) + a'^Zs - 1) V^S^(x) + a'^Zr V^T^(x) + a'^Z^ xi^)] 

+ E</2<f(^)- (19) 

Since it is not relevant to the present discussion we have left unspecified the exact lattice 
form of the SUSY current (point-split or local) and derivative (backward or symmetric). 
The same applies to the other dimension-9/2 operator appearing in Eq. (p!9|), namely 
the divergence of the mixing current T^{x). It may be defined in analogy to S'^(x) as 
point-split 

T!f'\x) = J]7.Tr(p^^j)(x)f/t(x)A(x + /i)f/,(x) + Plfj\x + f,)U,{x)X{x)Ul{x)) (20) 
or local 

T('°'=)(a;) = 25^7.Tr(p(:')(^)A(a;)) , (21) 

with the lattice derivative chosen as in the case of the current S^{x). From the above 
discussion obviously follows that different lattice currents S^{x) and T^{x) are associated 
with different values of Zs and Zj-. 

The last term on the r.h.s. of Eq. (|T9|) reflects the mixing of the operator Oii/2{x) with 
other bare operators O^^^^i^) of equal dimension. The reason Eq. ([19|) has been expressed 
in terms of the renormalized ones 0^^^^{x) will become clear shortly. The multiplicative 
renormalization Zii/2 and the mixing coefficients Z^y^ logarithmically divergent in 
perturbation theory. Solving Eq. ([T9|) for Ou/2{x) and substituting it in WI (0) one gets 

Zs ((V^5^(x)) Qiy)) + Zt {(VMx)) Q{y)) = ms (x(x)g(i/)) + 0(a) , (22) 

where the subtracted mass ms is given by 
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In deriving Eq. (^) we have relied on the vanishing in the continuum hmit of the corre- 
lation 

a ([Z~,),0^,,,{x) - J2 4/20!?/?(^)] Q{y)^ = 0{a) , (24) 

which is valid on-shell, x^y (recall that ^11/2, only logarithmically divergent). 

By using general renormalization group arguments (see e.g. 0) one can show that Zs, 
Zt and Z^, being power-subtraction coefficients, do not depend on the renormalization 
scale /i defining the renormalized operator (|T^. Consequently dimensional considerations 
imply Zs = Zs{go, t^oCl), Zt = Zt^Qo, rnoo), Z^ = Z^^qq, mQo). The requirement of a well 
defined chiral limit of the theory implies in particular that the dependence of Zs and Z^ 
on the gluino mass is vanishingly small in the continuum limit. In simulations at fixed 
lattice spacing this dependence is treated as an 0(a) effect. 

In QCD the lattice chiral WI, a relation analogous to Eq. (|2^) , leads to the definition 
of an axial current A^{x) = Z^A^If^x) where A^^*{x) is a generic discretization of the 
chiral current. A rigorous argument [Q, |^ shows that the current A^(a;) coincides with 
the correctly normalized continuum chiral current. It satisfies the appropriate current 
algebra. It is tempting to associate by analogy the quantity S^{x) = ZsSfj,{x) + ZtT^{x) 
with the correctly normalized continuum SUSY current. An attempt to reproduce in 
this case the QCD argument fails. This is because the proof in QCD relies on two key 
properties (cf. Sec. 3.1 of Ref. 0): 

1. The axial variation of the quark field is proportional to the field itself. 

2. The gauge fixing term is invariant under axial transformations. 

The above statements, valid for the axial symmetry, do not apply to SUSY. Thus, to 
the best of our knowledge, SUSY Wis cannot be used in a way analogous to the QCD 
chiral ones in order to prove the non-renormalization theorem for the current S^j,{x). 
Explicit one-loop calculations in lattice perturbation theory may shed some light on this 
issue [0, If the correctly normalized SUSY current coincides with S^{x) (or is related 
to it by multiplicative renormalization), it is conserved when ms vanishes. This is the 
restoration of SUSY in the continuum limit. 

4 Insertion operators 

In this Section we turn our attention to the insertion operator Q{x) of the WI (^). The 
operators V^S'^(x), V^T^(x) and x(x), which will be in the following collectively denoted 
as 'sink operators', transform according to the bispinorial representation of the Poincare 
group in the continuum. In order to get a non-trivial WI, the insertion operator Q{x) is 
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required to contain (at least) one non-zero spin-1/2 component. Thus, given a composite 
operator O which is a Majorana bispinor, Q{x) may be chosen to be of the form 

Q{x) = 0^{x) = C-^0{x) . (25) 

Clearly this operator must also be gauge invariant. 

We consider the zero spatial momentum WI obtained by summation over the spatial 
coordinates of Eq. 



J2 ((Vo5o(x)) 0^{y)) + ZrZ,' H^oUx)) 0^{y)) = 

X X 

X 

Note that in the above equation the three correlation functions are 4x4 matrices in 
Dirac indices. In numerical simulations these bare correlation functions can be computed 
at fixed lattice bare lattice parameters (3 = 2Nc/ Qq and n. Thus by choosing two elements 
of the 4x4 matrices, a system of two equations can be solved for rusZ^^ and Z-rZ^^. 

One must clearly ensure that these two equations are non-trivial and independent. 
To do this, let us consider the correlations containing the SUSY current (identical con- 
siderations apply for the other two correlations). Written explicitly with its Dirac indices 
this reads 

CiT^it) = E {(^oSoUx) Op{y)) , t = xo-y,. (27) 

X 

We consider dimensionless correlations, since these are the quantities actually computed 
in simulations. The above 4x4 matrix can be expanded in the basis of the 16 Dirac 
matrices F 

r 

Using discrete symmetries, see Appendix C, we can show that the only surviving contri- 
butions are (here Dirac indices are contracted) 

Ci^'^^Kt) ^ X]((^o(^)(^(2/))> (29) 

X 

C^rHt) = X]((V^o(^)7oO(2/))> . (30) 

X 

Due to the Majorana nature of the operators, the correlations ' (t) and C\f,' (t) 

are real. In conclusion for a given insertion operator we determine the dimensionless 
quantities amsZg^ and Z^Zg^ by solving the system of two equations 

Cf^)(t) + {ZtZs') C'^'''\t) = {amsZs')d^'''\t) 
C(f^)(t) + (Z2^Z^-i)C™(t) = iamsZs')c!yf''\t) . ^ ^ 



We now turn our attention to the choice of suitable insertion operators Q{x). Practi- 
cal considerations suggest the use of the lowest- dimensional insertion operators with the 
suitable symmetry properties. In our case this means the d = 7/2 gauge invariant bispinor 



Tr[p(:')(x)A(x 



(32) 



which is a tensor of 6 components in the Lorentz indices. Since the sink operators have 
spin-1/2, we must project out of the above Lorentz tensor the spin-1/2 components. 
Examples are So, Tq, x and 



(33) 



(only spatial plaquettes are taken into account). Since the Lorentz tensor of Eq. (0) has 
only two independent spin-1/2 components (see Appendix C for a detailed discussion), 
not all of the above operators can be independent. Indeed, they are related by 



7oTo(x) - 2x 



X] 



2^,{^Mx) - 2/'^\x)) 



(34) 
(35) 



We see that two independent systems of equations (3T) exist for two choices of dimension- 
7/2 insertion operators O. This redundancy can in principle be used in order to check 
lattice artifacts. 



5 Simulation of the model with light gluinos 

We simulate the N=l SU(2) SYM on a 12^ x 24 lattice at /3 = 2.3. This value of (3 
corresponds to the lower end of the approximate scaling region in pure SU(2) lattice 
gauge theory. In the full theory virtual loops of gluinos contribute to the Callan-Symanzik 
/5-function. The consequence is that, for fixed /3, the lattice spacing is decreased. 

The scaling properties of the model with dynamical gluinos were studied in detail in 



T0| . There values of k up to k = 0.1925 were considered. In that region of masses the 
observed effect coming from the dynamics of the gluinos was mainly the overall renormal- 
ization of the lattice spacing due to the fermionic virtual loops. The change of dimen- 
sionless ratios of masses and string tension where only moderate up to «;< 0.1925 where 
most of the simulations were performed. 



The set of configurations for the lightest gluino produced in |T0|, /t=0.1925, is taken 
here as a starting point. We further simulate the model at lighter gluinos, at k = 0.194 
and 0.1955. The largest of these hopping parameters coincides with the central value of 
the estimate of from the study of the finite volume gluino condensate, Kc = 0.1955(5) 
0. That determination was however obtained on a relatively small lattice (6'^ x 12) and 
the value of Kc is likely to be underestimated. In fact, anticipating results of the present 



study, the gluino mass starts decreasing significantly only for k > 0.194. For k = 0.1955 
the gluino mass is quite small but still appreciably different from zero. This is also evident 
in the simulation process. 

Simulating light fermions in a reasonably large physical volume is a challenging task 
from the algorithmic point of view. The difficulty is related to very small eigenvalues of 
the fermion matrix. The relevant parameter in this context is the condition number of 
Q^, where Q = 75(5 is the Hermitian fermion matrix^. For a given simulation volume 
the condition number gives an indication of the 'lightness' of the gluino. For the lightest 
gluino, at k = 0.1955, we had condition numbers O{10^). A direct comparison with the 
more familiar case of QCD is not possible since the simulation of SYM is generally less 
demanding. 

Another difficulty related to light fermions is the shrinking of the physical volume 
due to the renormalization of the lattice spacing by fermionic virtual loops. We expect 
that already at k=0.194 the low energy bound-states mass spectrum is strongly affected 
by the finite-size scaling on our 12'^ x 24 lattice. (See also Subsec. p.5| ). The situation is 
different for the main subject of this work, the SUSY Wis. The Wis hold also on a finite 
volume with volume-dependent coefficients. These are however essentially renormaliza- 
tions defined at the scale of the UV cutoff a~^. Our volumes should be consequently large 
enough for an accurate determination. 

All the numerical computations of this work were performed on the two 512-nodes 
CRAY-T3E machines at the John von Neumann Institute for Computing (NIC), Jiilich, 
with 307.2 and 614.4 GFLOPS peak-performance respectively. The CPU cost of the 
simulation was ~1 GFLOPS Year sustained (~3-10^^ f.p.o.) for each of the two simulation 
points of this work. 



5.1 The TSMB algorithm: simulation parameters 

The TSMB algorithm used in the simulations is defined in |jT2|, |l3l. The multi-bosonic 
updating with the scalar pseudofermion fields was performed by heatbath and overrelax- 
ation for the scalar fields and Metropolis sweeps for the gauge field; we refer to |]10[ for 
more details on the implementation of the TSMB algorithm in the case of the N=l SU(2) 
SYM. 

In Table |l] we report the parameters for the simulations performed in this work. We 
also include for reference the parameters for the simulation at k = 0.1925 performed in ||10|| . 
We briefiy explain the meaning of the various symbols (see also |TD|); 



[e. A] (columns 2 and 3) is the presumed domain of the eigenvalue spectrum of Q^; this 
is also the domain of validity of the polynomial approximation of the fermionic measure; 
^i,...,4 (columns 4 to 7) are the orders of the polynomial approximations used in the simu- 
^The condition number of a matrix is defined as the ratio between its largest and smallest eigenvalue. 
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Table 1: Parameters of the numerical simulations at (3 = 2.3. The run at «; = 0.1925 
was performed in |1^. The notation is explained in the text. 



K 


e- 10^ 


A 


n-i 

' " 1 




o 


n4 










updates 


offset 


' LCLl 


0.1925 


3.0 


3.7 


32 


150 


220 


400 


1 


3 


1 


1 


216000 


50 


9 


0.194 (a) 


0.8 


4.5 


38 


280 


320 


400 


1 


2 


8 


8 


2250 


50 


9 


0.194 (b) 


1.0 


4.5 


24 


160 


200 


400 


1 


2 


8 


8 


2700 


20 


9 


0.194 (c) 


1.0 


4.5 


24 


120 


160 


400 


1 


2 


8 


8 


1620 


20 


9 


0.194 (d) 


1.0 


4.5 


28 


160 


210 


400 


2 


2 


8 


8 


35460 


20 


9 


0.1955 (a) 


0.2 


5.0 


32 


420 


560 


840 


2 


10 


4 


4 


5040 


30 


8 


0.1955 (b) 


0.125 


5.0 


32 


480 


640 


960 


2 


10 


4 


4 


27672 


15 


8 


0.1955 (c) 


0.125 


5.0 


32 


480 


640 


960 


6 


6 


12 


12 


33120 


10 


8 



lation and measurement process; in particular, ni and n2 are the orders of the polynomial 
approximations in the local update and in the global accept-reject step (noisy correction) 
respectively; is the order of the polynomial used for the generation of the noisy vector 
in the noisy correction; n4 is the order of the polynomial used for the computation of 
the reweighting factors in the measurements; nHB,noB,nM (columns 8 to 10) indicate 
the heatbath, overrelaxation and Metropolis sweeps performed at each step of the local 
update; unc (column 11) is the number of Metropolis sweeps separating two consecutive 
global accept-reject steps; finally columns 12 to 14 report the total number of updates 
at equilibrium, the offset between measurements, and the number of independent lattices 
simulated. 



Table 2: Exponential and integrated autocorrelation of the plaquette measured in 
update cycles. In square brackets the integrated autocorrelation is measured in no. 
of Dirac matrix multiplications. The data at k = 0.1925 are taken from [^. 



K 


run 


^exp 


0.1925 




378(37) 


0.194 


d 


249(68) 


0.1955 


b 


220(50) 


0.1955 


c 


210(40) 


0.1955 


b,c 


260(30) 



675(200) [1.15(34) 10^] 

272(83) [0.75(23) 10^] 

280(70) [1.49(37) 10^] 

250(40) [1.71(27) 10^] 
420(50) 
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Figure 1: Distribution of the smallest eigenvalue of for the two simulation points of 
this work. The dashed line indicates the value of e used in the simulation. 



5.2 Autocorrelations 

Tuning the various parameters of the algorithm (e, A, ni,2,3, ^hb, nov, etc.) is essential to 
get an optimized updating. As an optimization criterion we have considered the autocor- 
relation of the plaquette. In particular, the order of the first polynomial ni was increased 
until the acceptance probability in the noisy correction reached ~50%. With this choice 
the two steps of the updating process equally contribute in shaping the distribution of 
gauge configurations. Larger values of ni have the effect of increasing autocorrelations 
with no substantial improvement of the algorithm. In the runs at /€ = 0.1955 the update 
of the pseudofermionic fields was performed by iterating twice a sub-sequence of nHs/'^ 
heatbath and uob/'^ overrelaxation sweeps. In Table ^ the integrated autocorrelations 



for the various runs are reported. The data for k = 0.1925 are taken from [T^. A better 
tuning of the parameters of the algorithm allowed to keep autocorrelations down at low 
levels in spite of an increasingly light gluino. 



5.3 Smallest eigenvalues and reweighting factors 

To monitor the accuracy of the polynomial approximation in the updating process, we 
constantly checked the smallest and largest eigenvalue of Q^. The distribution of the 
smallest eigenvalue for the two simulation points of this work is reported in Fig. |l|. The 
vertical dashed line indicates the value of e used in the simulation. 

Extremely small eigenvalues can exceptionally occur without substantial harm. The 
corresponding configuration would then be suppressed at the measurement level by the 
reweighting. We calculated reweighting factors for sub-samples of configurations. These 
turn out to be gaussian distributed with average ~ 1, see Fig. 0. For n = 0.194 we also 
observe a short tail towards small values. These distributions are consistent with the ab- 
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Figure 2: Distribution of the reweighting factors with gaussian fit. 

sence of extremely small eigenvalues in the ensembles. The effect of the reweighting turns 
out to be negligible compared to the statistical fluctuations for the quantities considered 
in this study. This confirms the overall accuracy of the simulation algorithm. 



5.4 The sign of the Pfaffian 



The fermionic measure implemented in the updating algorithm is given by \/ det((5). 
The actual measure for Majorana fermions is instead given by the Pfaffian of the anti- 
symmetric matrix M = CQ 



Pf(M) = Vdet(g) ■sign(Pf(M)) . 



(36) 



The previous formula implies that the configurations obtained by Monte Carlo updating 
should be further reweighted in the measurements by sign(Pf (M)) . This could potentially 
introduce difficulties. Indeed, were positive and negative signs almost equally distributed 
in our ensembles, a cancellation could occur in the statistical averages. The significance of 
the samples would then be close to zero. This occurrence is know as 'the sign problem'. In 
the N=l SYM with Wilson fermions the two signs are expected to be equally distributed 
for K > He- For k, < Kc, sign-flips should be suppressed if the volume is large enough (the 
situation right at the critical point is still unclear). 

We determined the sign of the Pfaffian for a subsample (10%) of configurations at 
fi;=0.194 on our 12^ x 24 lattice. We used the method explained in |ni| consisting in 
following the flow of the eigenvalues of Q as a function of the hopping parameter. A 
theorem ensures sign(Pf(M)) = 1 for small k's. The sign of the Pfaffian flips whenever 
an eigenvalue crosses zero in the flow. We never found such zero-level crossings, implying 
that sign(Pf (M)) is always positive for the considered sub-sample. Examples of eigenvalue 
flows are given in Fig. ^ We conclude that statistically less than 0.5% of the configurations 
have negative sign. Negative Pfaffians were detected in [IIO] on a smaller lattice (6^ x 12) 
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Figure 3: The flow of eigenvalues as a function of the hopping parameter for different 
configurations produced at j3 = 2.3 and k = 0.194. The vertical line indicates the hopping 
parameter of the simulation. 

for K > 0.196. Absence of sign- flip of the Pfaffian in our ensembles is also supported 
by the observation that extremely small eigenvalues do not occur and, consistently, the 
distribution of the reweighting factors does not extend to zero. Indeed each sign-flip 
under continuous modification of the gauge configuration would imply the crossing of a 
configuration with an exact fermionic zero-mode. 

5.5 Determination of rg/a 



In the scale parameter tq has been proposed as a reference scale for gauge configura- 
tions. It is defined by 

r^F(ro) = 1.65 (37) 

where F{r) is the force between static fermionic color sources in the fundamental repre- 
sentation. In our case the determination of ro/a allows us to monitor the reduction of the 
simulation volume with increasing k. 

We first calculated the static potential between color sources in the fundamental 
representation V{r); this is needed only at intermediate distances. The potential can be 
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Table 3: Determination of ro/a for the sets of configurations considered in this 
study. We also report the ratio L^jr^ where is the spatial lattice size. 





ro/a 


Lx/ro 


0.1925 


6.71(19) 


1.79(5) 


0.194 


7.37(30) 


1.63(7) 


0.1955 


7.98(48) 


1.50(9) 



estimated from Wilson loops W^(r,t) by 

V{v,t) = In 



{31 



Wir,t + 1] 

To reduce the noise we used APE smearing I^Bj on the spatial links with uape = 12 
iterations and e^PE = 0.5. Nevertheless the results get unstable for the larger Wilson 
loops, so we could not make an extrapolation to large t. We instead took the result for 
t = 2, calculated ro/a for these values and added the difference to the result from t = 3 
to the total error. 

To get ro/a from the potential we followed Their general procedure would 

be to make a fit to 



V{v) = Vo + ar - e 



+ f 



where 



An 



cos fk ■ r) 



(39) 



(40) 



removes lattice artifacts exactly to lowest order of perturbation theory. However it was 



not possible with our data to make fits to all four parameters. We therefore followed ||2J| 
and made a three parameter fit fixing e = 7r/12 and a two parameter fit fixing in addition 
/ = 0. One can now extract ro as 



1.65 -e 



(41) 



Our final estimate for ro comes from the three parameter fit, and the difference to the 
two parameter fit was again added to the Jackknife estimate for the error. 

In Table ^ we report the results for ro/a for the three sets of configurations considered 
in this study. From data in column 3 we see that the reduction of the physical lattice size 
is less than 20% for the lightest gluino in our simulations. 



6 Numerical analysis of the SUSY Wis 

In Section ^ we have shown how the quantities amsZ^^ and ZtZ^^ can be obtained 
by solving the system of equations (|3TD for a given insertion operator. We have argued 



15 



that there are two such independent dimension- 7/2 operators. We choose the operators 
x''"''\x) and To^'°'^(x) defined in (0) and (|21]) for our analysis. Discrepancies between 
results obtained with these two operators signal the presence of systematic effects due to 
the lattice discretization. 

In practice we use two methods to obtain amsZg^ and ZtZ^^. The first method is the 
most straightforward, consisting simply in solving the system of equations ( |3lD for each 
time-separation t. Results obtained in this way are reported in Figs. ^ and they should 
be independent of t when contact terms are absent for large enough time-separations. The 
second method, explained in detail in Appendix D, consists in constructing an overdeter- 
mined system of equations for several consecutive time-separations (tmm, ■ ■ ■ ,Lt/2) and 
fitting simultaneously for all time-separations. Results obtained in this way are reported 
in Tables |^-|^. Care should be taken in choosing t^j^ so as to avoid time-separations in 
which sink and insertion operators give rise to contact terms. 

Previous experience with the mass spectrum shows that point-like projection- 
operators give a poor signal for the correlations. This is mainly due to the gluonic content 
of both sink and insertion operators. The inconvenience is expected to be even more severe 
in the present case. Indeed sink operators involve time-derivatives which are subject 



to large statistical fiuctuations. The problem was solved in [|T0|, ^ by smearing the 
projection-operator for the gluino-glue bound state. Combined APE |^ and Jacobi ||2^ 
smearing were performed on the gluon and gluino fields respectively. In the present 
case we apply the same procedure only for the insertion operator. We searched for the 
optimal smearing parameters by analyzing sub-samples of gauge configurations. The 
set of parameters employed in the final analysis was: A^jacofej=18, -ft'jacofej= 0.2, Nape=^, 
eyiPE=0.5 (set A). In one case {k = 0.1925 and insertion operator we have also 

considered e^p^=0.1 and remaining parameters as in set A (set B). The APE- Jacobi 
procedure is not completely satisfactory in the case of insertion operators extended in 
the time-direction such as Tq°^\x). In our case with dynamical gluinos, a multi-hit 
procedure |2^ on the temporal links would not work since the noisy correction would be 
ineffective and anyway far too expensive. 

The inversions of the fermion matrix Q required for the computation of the cor- 
relations were performed by the conjugate-gradient method. The number of iterations 
necessary for a good accuracy on the final result increases for light gluinos. Convergence 
was improved by preconditioning the Hermitian fermion matrix. The residuum for the 
conjugate-gradient was chosen by requiring that the final accuracy on the determination 
of amsZg^ and ZxZg^ was <5 ■ 10~^. With this choice, ~1100 iterations were needed on 
average for the lightest gluino at k = 0.1955. The computing power used was correspond- 
ingly ~6 ■ 10^^ f.p.o.; this is about 20% of the amount employed for the generation of the 
gauge-fields. The site y of the insertion was chosen randomly for each configuration. We 
checked correlations in simulation time between propagators involved in the Wis. With 
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K = 0.1925 



Local Current 
K = 0.1925 



s 





Figure 4: anisZ^^ as a function of the time-separation t with insertion operator 

the random choice of y the correlations between two consecutively measured propagators 
turns out to be negligible (less than 0.05). Consequently a naive jackknife procedure can 
be used for the error analysis on amsZ^^ and ZtZ^^. 

At K = 0.1925 we have also used a version of the operator x^^^K^) defined as in 
Eq. (^) but with the lattice field tensor given by the simple plaquette 

= ^ iu,.i^) - KM ■ (42) 

The drawback of this definition is that the properties of transformation under P and T 
are not the same as in the continuum. However comparison of results obtained with the 
two definitions of x^'^^K^) gives an indication of the size of discretization errors. 



6.1 Results 

In Figs. ^ and ^ we report the determinations of arrisZ^^ and ZxZg^ respectively, as a 
function of the time-separation t. The left column refers to the point-split SUSY current 
Sjl'^\x) while the right one to the local current Sj!°^\x). The insertion operator is x^''^^K^)^ 
which generally gives the best signal. A behavior consistent with a plateau can be observed 
for time-separations t > 3 for which there is no contamination from contact terms. The 
signal is rapidly washed-out by the statistical fluctuations for large time-separations. 
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Figure 5: Z^Z^^ as a function of the time-separation t with insertion operator 

The insertion operator Tq"^\x) containing hnks in the time direction gives larger 
fluctuations than the time-slice operator This is probably related to the poor 

performance of the Jacobi-APE smearing on operators extended in the time direction. 

For a given insertion operator, ZtZ^^ is subject to larger statistical fluctuations than 
amsZg^ as can be seen from Fig. ^ (notice the different scale). We have no theoretical 
justification for this outcome. 

The results for amsZg^ and Zj^Z^^ from a global fit over a range of time-separations 
t > tmin are reported in Tables Data are obtained by solving the overdetermined 
linear system as explained in Appendix D. An equivalent procedure consists in performing 
a least mean square fit on arrisZg^it) and Z^Zg^it) for t > t^i^ taking time correlations 
into account. This second procedure gives results consistent with the first one, with 
X^/d.o.f. ~ 1 when contact terms are absent. For the final estimates we take tmin = 3 
for time-slice operators and tmin = 4 for operators extended in the time-direction. This 
choice ensures absence of contact terms. 

Discretization effects can be checked by comparing determinations obtained with the 
two independent insertions X^'P^x) and To^'°'^(x). For k = 0.1925 one can also compare 
between different definitions of (simple plaquette definition of the lattice field ten- 

sor, Eq. (|42|), and different smearing parameters). It should be recalled at this point that 
different discretizations of the currents give different values for Zs and Zt- Consequently 
data from different discretizations should not be confronted. 
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Figure 6: ZxZg^ as a function of l/zt witli tlie insertion operator x^^^^K^) (fiUed dia- 
monds) and Tq°^\x) (filled triangles). 

An interesting point is the dependence of ZtZ^^ and amsZg^ on the hopping pa- 
rameter K. This is reported in Fig. ^ and |^ respectively. Data refer to determinations 
obtained with tmin = 3 for insertion operator x^^^^K^) and tmin = 4 for insertion opera- 
tor Tq"'^\x). The definition of x^^'^H^) is the one with clover field tensor and smearing 
parameters of set A (cf. Tables ^-0). We see that the combination of renormalization 
factors ZxZg^ shows no appreciable dependence on k. We recall (see discussion in Sub- 
sec, that the latter is an 0{a) effect. Fitting these results (only for insertion x^^^K^)) 
with a constant in 1/k we obtain ZtZ^^ = —0.039(7), for the point-split current, and 
ZtZ^^ = 0.185(7) for the local current. The renormalization is surprisingly small in the 
case of the point-split current. An estimate of ZtZ^^ for the point-split current at /? = 2.3 
can be obtained from the 1-loop perturbative calculation in |]TB|]. At order one obtains 
ZtZ-' ^ = -0.074. 

In Fig. ^ the determination of anisZg^ is reported as a function of the inverse hopping 
parameter. The expectation is that anisZg^ vanishes linearly when n Kc- We see 
a clear decrease when k is increased towards Kc] arusZ^^ starts dropping abruptly at 
K > 0.194. We refer to data with insertion Given the weak dependence of the 

renormalization factors on the hopping parameter, the relative decrease of amsZg^ as a 
function of k should roughly compare with that of the gluino mass itself. In the case of the 
point-split current the latter almost halves when passing from k = 0.194 to fi; = 0.1955. In 
the case of the local current we get an even smaller mass at «; = 0.1955. Finally we obtain 
a determination Kc by performing an extrapolation to zero gluino mass from two largest 
K-values. The result is Kc = 0.19750(38) with the point-split current and Kc = 0.19647(27) 
with the local one. These values can be compared with the previous determination from 
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Figure 7: anisZg^ as a function of I/k with the insertion operator (diamonds) 
and T^''°^\x) (triangles). A linear extrapolation is also reported. The filled triangle 
indicates the determination of Kc from the first order phase transition of ||^. 



the phase transition Kc = 0.1955(5) 0. 



7 Summary and conclusions 

The present study shows that the extraction of the ratios anisZ^^ and Z^Zg-^ from the 
on-shell SUSY Ward identities is technically feasible with the computing resources at 
hand. The main technical difficulty (related to SUSY) is that high-dimensional operators 
with a mixed gluonic-fermionic composition must be considered, introducing relatively 
large statistical fiuctuations. This difficulty can be handled with an appropriate smear- 
ing procedure. The non-perturbative determination of the ratio amsZg^ can be used for 
a determination of the critical hopping parameter Hc corresponding to massless gluinos. 
This can be compared to the independent determination of Hc, based on the chiral phase 
transition [Q. The fact that the two determinations are not in full agreement can be 
explained by the presence of systematic effects. These are predominantly 0(a) effects 
in our case and finite volume effects in the determination of Ref. The discretization 
error in the present approach can be checked by comparing results from the two inde- 
pendent insertions. The statistical uncertainty is however relatively large (10-40%) for 
the insertion operator involving links in the time-direction. Our results for arris Zg^ and 
Z^Zg^ are consistent with the Wis (|2^) with 0(a) effects comparable to the statistical er- 
rors. The non-perturbative results here presented will be complemented by an analytical 
perturbative calculation [|T^. 

The overall conclusion of our study is that lattice SUSY Wis can be implemented in 
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the non-perturbative determination of the correctly subtracted SUSY current and provide 
a practicable method for the verification of SUSY restoration in this framework, once the 
proper SUSY current is identified and the usual lattice artefacts (most notably finite cutoff 
effects) are kept under control. 
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Table 4: Summary of the results for amsZg aX (3 = 2.3 with point-split currents. 





operator 


'^min 




'^min ^ 


0.1925 




0.176(5) 


0.166(10) 


0.135(14) 


n 1 no c 
(J. 1925 


, ,(sv\ f 


0.1o2(dj 


0.152(11) 


0.150(16) 


0.1925 




0.1969(47) 


0.168(9) 


0.136(14) 


0.1925 


-'o 




0.132(16) 


0.124(21) 


0.194 




0.148(6) 


0.130(11) 


0.146(21) 


0.194 






0.095(27) 


0.090(27) 


0.1955 




0.0839(35) 


0.0820(7) 


0.053(14) 



* With plaquette field tensor. 



** With plaquette field tensor and smearing parameters B. 



Table 5: Summary of the results for amsZg^ at /5 = 2.3 with local currents. 





operator 




tmin — 4 


^min 5 


0.1925 


■^{sp) 


0.166(6) 


0.166(11) 


0.146(16) 


0.1925 


^(sp) (*) 


0.173(6) 


0.155(11) 


0.135(19) 


0.1925 


(**) 


0.1821(47) 


0.173(11) 


0.154(19) 


0.1925 


rp{loc) 




0.144(18) 


0.143(25) 


0.194 


^(sp) 


0.124(6) 


0.126(12) 


0.142(24) 


0.194 


rp{loc) 




0.076(30) 


0.098(35) 


0.1955 




0.0532(40) 


0.064(8) 


0.047(15) 



* With plaquette field tensor. 

** With plaquette field tensor and smearing parameters B. 
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Table 6: Summary of the results for ZtZ^ at /? = 2.3 with point-split currents. 



K 


operator 


tmin 3 


tmin 4 




0.1925 




-0.015(19) 


-0.036(31) 


0.045(56) 


0.1925 


^(sp) (^*^ 


—0.044(16) 


—0.096(33) 


0.01(6) 


0.1925 




-0.058(14) 


-0.044(32) 


-0.07(5) 


0.1925 


rpiloc) 




0.11(7) 


-0.03(7) 


0.194 




-0.038(19) 


-0.024(43) 


-0.08(7) 


0.194 


rp{loc) 

-'o 




0.11(13) 


0.02(13) 


0.1955 


■^{sp) 


-0.051(13) 


-0.064(26) 


-0.05(5) 


* With plaquette field tensor. 



** With plaquette field tensor and smearing parameters B. 



Table 7: Summary of the results for ZtZ^^ aX (3 = 2.3 with local currents. 



AC 


operator 


tmin 3 


tmin 4 


^min 5 


0.1925 




0.183(14) 


0.207(27) 


0.19(5) 


0.1925 


^isp) (*) 


0.176(14) 


0.184(28) 


0.21(5) 


0.1925 




0.146(11) 


0.159(25) 


0.139(45) 


0.1925 


q-l{loc) 




0.29(6) 


0.22(6) 


0.194 




0.202(15) 


0.176(33) 


0.186(6) 


0.194 


rp{loc) 

-'o 




0.27(9) 


0.30(11) 


0.1955 


-^{sp) 


0.179(10) 


0.170(21) 


0.170(45) 



* With plaquette field tensor. 

** With plaquette field tensor and smearing parameters B. 
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Appendix 



A Discrete symmetries 

The following Hermitian representation of the Euclidean 7-matrices is adopted: 
with anti-commutation property: 

{7;„7.} = 2V. (44) 

The matrix 75 is defined as 

/ 1 , 

75 = 7i72737o = \ ^ ^ ] , (45) 



-1 



and the anti-Hermitian matrix cr^j, reads 



(T^,u = ^b^„lu] ■ (46) 

A.l Parity P 

Denoting with xp the transformed coordinates, xp = (xq, — xi, —X2, —X3), the field trans- 
formations are: 

AP(x) = 7oA(xp) (47) 
X^x) = Xixpho (48) 
U^{x) = Uoixp) (49) 
U[ix) = uKxp-l). (50) 

A. 2 Time reversal T 

In this case xj = (— xq, xi, X2, xs), and the field transformations are: 

A^(x) = 7o75A(xt) (51) 
A^(x) = A(xt)757o (52) 
Ujix) = Uiixj-O) (53) 
Ujix) = Uiixj) . (54) 

(55) 

The clover-symmetrized field tensor P^^'^ (x) of Eq. ([T3|) has the same transformation 
properties under P and T as its continuum counterpart F^^{x). This does not apply to 
the simple plaquette lattice field tensor Pliu\x) of Eq. (^). 
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A. 3 Charge conjugation C 

The invariance of the gluino field under charge conjugation C (Majorana condition) reads 

A(x) = A^(x) = Cfiyx) , 

A(x) = A^(x) = X^{x)C-^ . (56) 
The spinorial matrix C is defined as 

C„,.= f 7 ° ) . (57) 

with the properties 

= -C = C^ (58) 

and 

C-'^,C = -7j , C-'cr,,C = -al , C-^C = f . (59) 



B Renormalization of Oii/2(x) 



We discuss here the renormahzation of the composite operator Oii/2{x) defining the con- 
tinuum hmit of X{x) according to power-counting P], see Eq. ([T9|) . In this context the 
particular lattice form of the operators is immaterial and continuum notation will be used 
for simplicity. 

In the present work we consider gauge invariant correlation functions of the opera- 
tor Oiii2{x). In this case one can restrict the analysis to mixing with gauge invariant 
operators. These operators must have in addition the same transformation properties 
as Oii/2(a;) under the hypercubic group and the discrete symmetries C, P and T. Fi- 
nally, dimensional considerations restrict the search to operators of dimension d < 11/2 



(dimension- 11/2 or higher mixing can be neglected as clarified in Sec. 3.1 
The operator Oii/2{x) transforms under 0(4) like the gluino field 

K{x) '^Saf3{R)\p{x) S{R) = exp{-i ^cj^^o-^^} . (60) 

The connection with 0(4) is given by 

S{R)-\S{R) = J2R>^u7u . (61) 

The coefficient |JJ^J^v represents the rotation angle in the plane (/i, v)\ the hypercubic group, 
a discrete subgroup of 0(4), is obtained by restricting these angles to multiples of 7r/2. 

Gauge invariant operators are obtained by taking the trace in color space of products 
of fields in the adjoint representation. Available operators of this kind are the gluino field 
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itself, the field tensor F^u{x) and the covariant derivative in the adjoint representation 
D^. Dirac and tensorial indices should be combined in such a way that the resulting 
operator transforms like (0) under hypercubic transformations. By using the algebra of 
the F-matrices 

^ = {]l,7^7^,7^7^,a^J (62) 
pairs of gluino fields can be assembled into bilinear expressions of the form 

0^^'^"(x) = A(x)r... \{x) (63) 

where Dirac indices are contracted and the dots indicate possible tensorial indices of the 
F-matrix. The transformation properties of the bilinears ( p^ ) depend on the F-matrix. 
They belong to irreducible representations of both 0(4) and the hypercubic group, see 
e.g. [H, 



F = 1,7^: scalar, /(+'") {ri^}) , 

r = 7M,7'7M: vectorial, (1/2, l/2)(+'-) (r^J; 

F = a^^: tensorial, (1,0) © (0, 1) (rf^) . 



The representations of 0(4) are indicated in the notation of [^]. In parenthesis we report 



the corresponding representations of the hypercubic group in the notation used in ||29 
(see the following). 

General classification. The general operator can be classified according to the 
number of contained gluino fields A(x), field tensors Ffj^i^{x) and covariant derivatives D^, 
n\, Up and n^). In the case of an even number of gluino fields, one can exploit the 
algebra of the F-matrices and build bilinears of the form ( |63D with no unpaired gluino 
field. Consequently all Dirac indices are contracted. Products made up by these bilinears 
0^*'*"'(a;), field tensors Ffj_^{x) and covariant derivatives transform according to a gen- 
eral tensorial representation of the hypercubic group O^^^...^^^. An example is given by 
the dimension-4 operator 

O^,^, =Tr[D^,(A(x)7^,A(a:))] . (64) 

The decomposition into irreducible representations of the hypercubic group of tenso- 
rial representations has been studied in Of course the representation ( |60D cannot 
be contained in these representations.^ One can consequently restrict the investigation 
to operators containing an odd number of gluino fields. An additional restriction to the 
possible structures of mixing operators comes from the general relation 

Tr [D,, ■ ■ ■ D,^Aix)] = d,, ■ ■ ■ 9^„Tr [A(x)] , (65) 



■^This is evident for example if one considers that a rotation of 27r according to (^0|) is not the identical 
transformation as in the tensorial representations. Relation ( |60|) defines a double-valued irreducible 
representation of 0(4) and of the hypercubic group. 
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holding for a generic adjoint field A{x). Since Tr[A(2;)] = relation (|65D excludes cases 
where nx = 1 and np = 0. This leaves only three possibilities for d < 9/2: 

a) nx = 3, Up = , riD = {d = 9/2) , 

h) nx = I, Up = 1 , riD = 1 {d = 9/2) , 

c) nx = I, Up = 1 , ud = {d = 7/2) . 

B.l nx = 3, np = , no = 

The most general gauge invariant three-gluino operator can be expressed as 

0«,...(x) = Tr [(A(x)r . A(x)) (F A(x)) J , (66) 

where as usual the dots indicate possible tensorial indices. Dirac indices are contracted in 
the bilinear. The transformation rule in the free Dirac index a is the right one given by 
Eq. (|60|). This Dirac index will be suppressed in the future notation. Oa,..Xx) defines a 
tensorial representation of 0(4) (and of the hypercubic group) in the remaining tensorial 
indices. Imposing invariance under 0(4) is equivalent to requiring that these indices are 
contracted so as to obtain a scalar. The procedure is standard and the result is 



Os{x) 
Op{x) 

Ov{x) 
Oa{x) 
Ot{x) 



Tr [(A(a;)A(a;))A(x)] 

Tr [(A(x)7^A(x))7^A(x)] 



Tr 



Tr 



Tr 



^(A(x)7/,A(x))7/,A(x) 

. 

^(A(x)7S,A(x))7^^A( 

. 

^(A(x)(T^^A(x))o-/,^A(x) 



(67) 
(68) 

(69) 
(70) 
(71) 



For Majorana fermions any three-gluino operator complying with 0(4) invariance, (|67|) - 
(|7l]), vanishes. The Majorana condition (^) can be directly used to show vanishing of 
Os{x), Op{x) and Oa{x): 



Fierz rearrangements 



Os{x) = Op{x) = OAix) = . 



Ov{x) = -Oa{x) = Os{x) - Op{x) 
Ot{x) = -Os{x) - Op{x) 



(72) 



(73) 
(74) 



imply vanishing of the remaining two operators Ov{x) and Ot{x). 
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The discussion is however not complete since the true symmetry of the lattice is the 
hypercubic one. So the question arises whether tensorial indices can be combined in ( |66D 
in a different way from (|69|)-([7TD, while still complying with the hypercubic invariance. 
The resulting operator would represent a potential Lorentz-breaking term in the renormal- 
ization of Oii/2{x). The argument can be made more rigorous by considering that Eq. ( pBD 
defines tensorial representations of the hypercubic group. The question is whether singlet 
representations of the hypercubic group, which are not 0(4) scalars, are contained in the 
tensorial representations (|55|). We rely in this on the detailed discussion of the subject 
contained in ^9 . 



Singlet representations of the hypercubic group which are not necessarily 0(4) scalars 
are contained in tensorial representations with even number of indices n > 4. In the 
classification ([66| ) the only possible candidate is the operator with T^^, = T'^^ = a^i, 
containing four tensorial indices 

(x) = Tr [{X{x)af,^^2X{x))af,.^^^X{x)] . (75) 

This means that 0(4)-breaking versions of the operator Ot{x), Eq. (^), are in principle 
possible. In order to decide the question we consider the general case of an operator with 
four tensorial indices Ol^^f_i^^.^^4^{x). One-dimensional representations of the hypercubic 
group which do not coincide with 0(4) scalars are given by p 



OfoU}!^) ri'V (77) 



The symbol {...} in ([77|) defines complete symmetrization on the Lorentz indices. In the 
case we are interested in, (0^^)^i^2M3Ai4(-^)' antisymmetry of a^,^ implies trivially vanishing 
of these new combinations of indices (combination (1771) would have the wrong symmetry 
anyway). Consequently no three-gluino operators complying with hypercubic invariance 
can be built. 

B.2 nx = 1, Up = I , nj) = 1 

Gauge invariance imposes the form 

0(x)^„^„^3,... = Tr [D,,{F,,,,{x)r,„ X{x)}] . (78) 

The symbolic expression D^{ABC ■ ■ ■} is a collective representation of all the ways the 
covariant derivative can act on the products of the adjoint fields A,B^C . . .. Since 
the adjoint covariant derivative satisfies the Leibniz rule, it is sufficient to consider only 
operators where acts on just one field A, B,C . . .. Again, tensorial indices must be 
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combined as to obtain singlets under the hypercubic group. 0(4)-invariance and P leave 
only two possibilities 

Tr[D4F^,(x)7.A(x)}] (t) (79) 

TT[e^i,paD^{F^p{x)j5ja\{x)}] {ii) . (80) 

Also in this case new hypercubic-invariant combinations could come only from 4-index 
tensors, e.g. 

C(a;)^i,^2,;.3,;.4 = Tr [L>^, {F^^^g (a;)7^, A (x) }] . (81) 

An analysis analogous to the one performed in the previous Subsection leads to the con- 
clusion that no new, non-trivial operators arise from combinations (f^) , ([77D . In this case 
antisymmetry of F^y{x) plays the key role. So we concentrate on the Lorentz-conserving 
operators ([79| ) and (pO|). 

Case (i). The first possibility according to the Leibniz rule for the covariant derivative 
in (ED 

Tr[(D^F^,(x))7.A(x)] (82) 

vanishes on-shell. Indeed, using the equation of motion for F^y{x) it can be rewritten as 
a three-gluino operator, vanishing identically according to the previous discussion. The 
second possibility for the Leibniz rule 

Tr [F^,{x)-i,D^\{x)] (83) 

is equivalent on shell to 

Tr[D,(F^,(x)7.A(x))] , (84) 

given the vanishing of (P2|). Using relation ( |55D we arrive at the operator 9^T^(x). 
Case (ii). This goes along the same lines as case (i). The relation 

^^lupcTl^l<T = -crupl^l + {^^plu - ^^iulp) (85) 

can be used. The part containing the 7^ matrices reduces to the case (i) already consid- 
ered. The new combination is: 

1i[D^{F,p{x)a,p^^\{x)}] . (86) 
Again, there are essentially two possibilities: 

Tr[F^,(x)a^,^A(x)] (87) 

and 

i:i[D^{F,p{x)a,p-i^\{x))] . (88) 

On-shell, the first operator reduces to the lower dimensional operator x{^) (Q) 
equation of motion for the field A(x), and can be neglected in this discussion. Rule (|65D 
shows that the second operator is just the divergence of the SUSY current dp,Sp{x). 
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In summary, the systematic scan of all possible dimension- 9/2 operators with the 
required symmetry properties and the on-shell restriction leaves us with the two operators 
d^Sf,{x) and d^T^{x). 

B.3 n\ = 1^ np = ^ , no = 

Using gauge invariance one gets of the general form 

0(x)^„^,,... = Tr[F^,^,(x)r.A(x)] . (89) 

The only hypercubic-invariant combination of indices is 

TT[F^,{x)a^,X{x)] = x{x) . (90) 

Again, no Lorentz-breaking combination appears. 

This exhausts all possibilities for the power subtractions of Oii/2{x) and implies the 
form (ID of the SUSY Wis. 



C Insertion operators 

In this Appendix we use for simplicity, as in the previous one, notions of the continuum. 
Even if the true symmetry of the lattice is the hypercubic one, the analysis is carried over 
by using the more restrictive Lorentz invariance. Indeed Lorentz breaking terms should 
not be considered, given that 0(a) effects are neglected. 

Consider an insertion operator transforming under a generic representation of the 
Lorentz group (we restrict the discussion to operators depending on one coordinate). Due 
to the spinoral character of the SUSY Wis we need spinorial operators. Operators of this 
type are 

0^,...f,„,a{x) = r^i...^„(x)^„(x) , (91) 

where T^i...^„(x) is an operator transforming according to the tensorial representation of 
the Lorentz group and ipaix) is a bispinor. Irreducible representations are projected-out 
by suitable (anti)symmetrizations, extraction of traces on the free indices. We consider 
zero spatial momentum Wis (p6[). These are obtained by taking the vacuum expectation 
value of operators like e.g. 

Of,^..,^„^af3{xo,yo) = Y j dxdy{doSo)aix,xo)Of,^,„f,„^f3{y,yo)- (92) 

The total angular momentum J of the above operator results from the composition of 
the spins of the two operators / dx {doSo)a{x,Xo) and / rfx 0^j...^„^(x, ?/o); the 'orbital' 
angular momentum is zero because of the summation on space coordinates. In order to 
get non-trivial Wis one must have of course J = 0. Since the operator d^Sij^{x) is a 
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bispinor, the condition J = implies that the insertion operator must contain at least 
one spin-1/2 component. The irreducible representations of the Lorentz groupQ with this 
property are of the type {{S, S') © {S', S)), with |^ - S'\ = 1/2. 

The lowest dimensional gauge invariant operator has dimension d = 7/2. It has 
the form ( PT| ) with tensorial part given by and ipi^) = A(x). Its transformation 

properties under the Lorentz group are given by 

((1, 0)0(0, 1))®((^, 0)0(0, ^)) = ((^, 0)0(0, ^))©(^, ^)©(^, ^)©((i, 1)©(1, ^)) . (93) 

We see that two (and only two) representations on the r.h.s. of the above decomposition 
contain spin-1/2, namely ((|, 0)©(0, ^)) and ((|, 1)©(1, i)). The first component is given 
by the operator x{^)y the second is present for example in S'^(x) and T^(a;) transforming 
hke 

^) ® ((^, 0) © (0, i)) = ((i, 0) © (0, i)) © ((^, 1) © (1, i)) . (94) 

Pure spin-1/2 operators are given by x(^)) So{x) and To{x), while Si{x) and Tj(x) con- 
tain also spin-3/2 components. Eq. (0) implies that any spin-1/2 component of any 
dimension-7/2 operator can be expressed as a linear combination of two operators chosen 
as a basis, for example x{^) Tq{x) (see e.g. relations (|51D, (|5BD). Consequently only 
two operators give independent Wis at the lowest dimension. 

For a given insertion operator the Wis consist of a number of independent equations 
which equals the number of rotational invariant (spin-0) components contained in the 
operatorial product (^2]). A straightforward analysis including also P reveals that these 
are two for O = x(x), S'o(x), To(x), as given in (p9D and (|30D . It should be stressed 
that considering other insertion operators like for example Si{x) or Ti{x) would not give 
additional information to that given by e.g. x{x)y To{x). One would get only different 
combinations of the same on- shell Wis. 

The discrete symmetries T and C imply relations for the correlations involved in the 
SUSY Wis. For the correlations defined in (|2^), (l30|) these are 



T: C(^-^)(t) = 7o75C(^'^)(-t)757o (95) 
C: C(^'^)(t) = 7o7275(C(^'^))*(t)75727o • (96) 
In terms of the components defined in (^), relation (|95|) reads 

Ct°>(-t) = Cf 

C<f°)(-t) = -C<f°)(i) . (97) 

In the case O = So,Tq an extra minus sign must be included. Relation (|96|) ensures reality 
of the two components 

Cf°)(t) = (Cf'^^(t))* 
ag;^(t) = (C(J'^)(t))* . (98) 

*As usual we take into account also P. 
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Properties related to P and T apply unchanged for the lattice theory if the clover- 
symmetrized lattice field tensor pIw\x) is used for the insertion operators. The simple 
discretization Pj^y\x) breaks P and T and additional 0(a) pseudoscalar components are 
present in the correlations. In this case, for operators extended in the time-direction, the 
time-reflection properties ( P7D are violated. 



D Linear fit 

We perform a linear fit to solve the WI (|3TD for ZxZ~^^ and amsZ~^^ including several 
consecutive time-separations (tmm, ■ ■ ■ ,Lt/2). We define 

A = ZtZs\ B = amsZs^ , (99) 

and Xi^t, yi,t and zi^t as the different components of the correlation functions at different 
times (see Eqs. (PP|),(PD|)) 



X2,t — (--70 l^J ? y2,t — ^-yo ; Z2,t — <-^70 l^J 



The overdetermined system reads 

Xi^t + Ayi^t = Bzi^t, i = l,2, t = tmin, ■ ■ ■ , Lt/^ . (101) 
We get the best estimates for ZtZ^^ and amsZ^^ (A and B) by minimizing the quantity 

2 Lt/2 
^—1 i — imin 



With the conditions 

dH 
'dA 

dH 
'dB 

and the definition 
A and B are given by 



2 J2 i^i,t + Ayi,t - Bz,^t) = (103) 

i,t 

-2 ^i.t {^i,t + Ayi,t - Bzi^t) = (104) 



i,t 



X] Htyi,t =<x,y>, (105) 



i,t 



^ ^ <y,z >< x,z > - < x,y >< z,z > 

<y,y >< z,z > - <y,z 

< X, z >< y,y > — < x,y >< y, z > 

B = — . 107 

< y,y >< z,z > - < y,z >^ 
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